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Two- and three-dimensional electron gases with a uniform neutralizing background are studied at negative 
compressibility. Parametrized expressions for the dielectric function are used to access this strong-coupling 
regime, where the screened Coulomb potential becomes overall attractive for like charges. Closely examining 
these expressions reveals that the ground state with a periodic modulation of the charge density, albeit expo- 
nentially damped, replaces the homogeneous one at positive compressibility. The wavevector characterizing 
the new ground state depends on the density and is complex, having a positive imaginary part, as does the 
homogeneous ground state, and real part, as does the genuine charge density wave. 



I. INTRODUCTION 

The nature of the metallic state of a dilute two-dimensional 
electron gas (2DEG) in high-mobility silicon MOSFETs, first 
observed by Kravchenko et aln nearly a decade ago, has 
still not been established. Various scenarios have been pro- 
posed to explain the new conducting state, ranging from su- 
perconductingj-to non-Fermi liquid, to percolation, to classi- 
cal (see Refs.Btl for overviews and references to the original 
literature). The unique character of this state is reflected by 
its thermodynamic signatures of both negative pressure and 
compressibility!! 

Although recognized as a theoretical possibility for quite 
some time — at least in three dimensions!!], an electron gas with 
a uniform neutralizing background was expected to be unsta- 
ble at negative compressibility, and the model meaningless. 
Later theoreticalQ and also numerical!] studies nevertheless 
considered both the region with positive and negative com- 
pressibility without encountering conceptual difficulties. 

Eisenstein, Pfeiffer, and West§ more recently concluded 
that the global features of their 2DEG compressibility data 
to very low densities, where the system appears to undergo a 
quantum phase transition to an insulating statetj, could be ex- 
plained simply using the Hartree-Fock approximation, which 
consists of including the exchange contribution to the ground- 
state energy of a free fermion gas. Their analysis does not 
account for impurities, which are of paramount importance at 
and beyond the metal-insulating (MI) transition. 

Hani et a/.i reached similar conclusions based on their com- 
pressibility data of a dilute two-dimensional hole gas (2DHG). 
By placing several single electron transistors directly above 
the 2DHG, they could in addition determine possible spatial 
variations in the system. They concluded that the metallic 
state at negative compressibility is homogeneous in space and 
well described by the Hartree-Fock approximation. As the 
system crosses to the insulating state and the Hartree-Fock 
approximation ceases to be applicable, they found the system 
to become spatially inhomogeneous. They interpreted this ob- 
servation as supporting scenarios in_wluch the MI transition is 
described as a percolation processtMCJ 

Using density functional theory in the local density approx- 
imation, Shi and XieLLa recently investigated the spatial dis- 
tribution of the electron number density of a 2DEG in the 
presence of impurities. As usual, impurities were included by 



coupling the particle number density to a fluctuating potential 
with a random distribution. The parameter characterizing the 
Gaussian distribution, which we will refer to as the impurity 
strength, determines the roughness of the impurity landscape. 
Their numerical study incorporates the Monte Carlo data on 
the ground-state energy of a clean system by Tanatar and 
CeperlyQ. For given impurity strength, the study shows that at 
some average density n islands of very low densities form in a 
metallic sea of high densities. At higher average densities, the 
sea level is high enough to fill all of the valleys of the impurity 
landscape, and the system is homogeneous. With decreasing 
n, the sea level drops and the insulating islands grow. At a 
critical value, the islands percolate the system, which at this 
percolation threshold becomes insulating. The sea level now 
dropped to the extent that the electrons are confined to the val- 
leys of the impurity landscape. According to this picture, the 
transition to the insulating state takes place at lower densities 
for cleaner systems. This is in accord with experiments, where 
the lowest critical density, n c = 7.7±0.4x 10 9 cm~ 2 , was ob- 
served in an exceptionally clean 2DHSE1 Shi and XieEl also 
studied the importance of the Coulomb interaction by com- 
paring with the case where this interaction is turned off. They 
found that at an average density where the interacting system 
was still metallic, the free electron gas, which is insulating, 
had formed a few isolated lakes of high density at valleys of 
the impurity landscape. 

In this paper, we postulate on the nature of the ground state 
of electron gases at negative compressibility. In doing so, we 
ignore impurities, as justified by the experimental findings that 
this metallic state is homogeneous in space and that the com- 
pressibility data is already well described by the Hartree-Fock 
approximation of a clean electron gas. The study emphasizes 
the three-dimensional electron gas (3DEG) as much more in- 
formation from theoretical studies is available than for the 
2DEG. Most of the 3D conclusions, however, also apply to the 
2D case discussed in Sec. VI. The following two sections fo- 
cus on the changes brought about when going from positive to 
negative compressibility. One of the most surprising changes 
is that a test particle acquires a screening cloud in its imme- 
diate vicinity which overcompensates the test charge, so that 
the Coulomb interaction becomes attractive for like charges. 
The changes are also discussed from the perspective of Fermi 
liquid theory, describing the electron gas after the screening 
mechanism has taken effect and resulted in only short-range 
interactions. The main result of this paper is contained in Sec. 
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[Tvj where it is argued that the ground state of an electron gas 
at negative compressibility is a charge density wave (CDW), 
although exponentially damped. The phenomena considered 
here are by no means specific to an electron gas, but also ap- 
pear in a charged Bose gas as discussed in Sec. [v|. 

II. SCREENING VS. OVERSCREENING 

A 3DEG is characterized by the screening of charges. 
When a test particle of charge Z is placed at the origin, 
p(x) = Z6(x), the system responds by rearranging its charge 
distribution to screen the external charge. The charge density 
Pind thus induced, 
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is determined by the dielectric function e(q, 0) at zero fre- 
quency, which encodes the static screening effects of the 
3DEG. The total induced charge 



d 3 Xp in(i (x) = Pind(q = 0) 



(2) 



cancels the charge of the test particle since the screening fac- 
tor l/e(q, 0) vanishes at zero wavevector q = |q|. The test 
charge is therefore always perfectly screened. However, the 
induced charge density differs in how it is distributed, de- 
pending on the electron number density n, or, equivalently, 
the ratio r s = a/ao of the average interparticle distance 
a = (3/47rn) 1 / 3 to the Bohr radius ciq = fi 2 /me 2 . (For va- 
lence electrons in 3D metals, r s , characterizing the strength of 
the Coulomb interaction, ranges fromS 1.8 to 5.6.) 

At weak coupling (r s < 1), corresponding to high elec- 
tron number densities, a screening cloud of opposite charge 
surrounds the test charge, and the screened Coulomb poten- 
tial decreases exponentially with increasing distance. On its 
tail, far away from the test charge, the potential has superim- 
posed a small oscillatory modulation of wavevector q = 2kp, 
with fcp the Fermi wavevector. These Friedel oscillations, 
leading to a periodic change in sign of the Coulomb poten- 
tial, originate from a singularity in the dielectric function at 
q = 2fcp, where electron-hole excitations start to develop an 
energy gapQ. 

A 3DEG's response to a test charge fundamentally changes 
at larger values of the coupling constant r s because of a quali- 
tative change in the dielectric function. Namely, at some value 
r s = f s , the-jdielectric function becomes negative for small 
wavevectorslii Rather than surrounded by a screening cloud 
of opposite charge in its immediate vicinity, the test charge 
now becomes overscreened. The Coulomb potential rapidly 
drops below zero with increasing distance and becomes at- 
tractive for like charges. Further away from the test charge, 
the potential exhibits an exponentially damped oscillatory be- 
havior and periodically changes sign similar to Friedel oscil- 
lations. The initial drop of the Coulomb potential to negative 
values and the resulting overscreening is, however, unrelated 
to Friedel oscillations^. Overscreening in a plasma was first 
noted in a 2D Bose gas with a 1 jx Coulomb potential^. 



The fundamental change in the screening behavior of a 
3DEG is paralleled by a modification of its ground state. At 
the critical electron number density, where the dielectric func- 
tion becomes negative for small wavevectors and a test charge 
overscreened, the compressibility changes sign as well. As 
argued in the following, the homogeneous ground state of 
the regime with positive compressibility then gives way to a 
ground state with a periodic modulation of the charge density, 
which is, however, exponentially damped. 



III. NEGATIVE COMPRESSIBILITY 

At weak coupling, .the dielectric function can be written for 
small wavevectors asQ 



lime(g,0) = l-u(g)xsc(0,0), 

g^O 



(3) 



with v(q) = 4ire 2 /q 2 the Fourier transform of the (un- 
screened) Coulomb potential, and x sc (0,0) the screened 
density-density response function at zero wavevector and fre- 
quency. The screened response function, a theoretic con- 
struct to be distinguished from the physical one, measures 
the response to a screened external field rather that the ex- 
ternal field itself. Physically, this function describes a ficti- 
tious system with only short-range interactions that are rem- 
nants of the long-range Coulomb interaction after the screen- 
ing mechanism of the 3DEG has taken effect!]. Derived from 
the Coulomb interaction, the short-range interactions of the 
fictitious system vanish in the limit r s — > 0. 

In writing Eq. (^]), the 0(q 2 ) term in the Taylor expansion 
of Xsc 0) is assumed to be substantially smaller than 1. This 
is true only at weak coupling. For example, when r s = 1, the 
value of the constant term in Eq. (|^) is reduced already by 
about 6 %. 

The compressibility sum rule relates the screened response 
function to the compressibility k of the 3DEG viaO 



l- / n\ ® n 2 

limY 8c (g,0) = -— = -n k, 

q->0 dp 



(4) 



with p the chemical potential. If positive, the compressibil- 
ity can be expressed in terms of the speed of sound c in the 
fictitious system with only short-range interactions as 



riK = 1/mc 2 . 
Equation (||) can then be cast in the equivalent form 



pl 



lim e{q, 0) = 1 , 
g^O c z q z 

corresponding to the spectrum oj 2 = lo 2 



pi 



(5) 



(6) 



c 2 q 2 of the 



plasma mode, where tt> p i denotes the plasma frequency, ui^ — 
Aime 2 /m. The plasmon spectrum differs from the gapless 
spectrum lo 2 = c 2 q 2 of the sound mode of the fictitious sys- 
tem in that, due to the long range of the Coulomb interac- 
tion, the former has an energy gap. At short wavelengths 
(cq > w p i), the difference is negligible, thereby allowing use 
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of the plasmon to access, via c, the screened interaction as a 
function of the coupling constant — at least for r s < 1. 

In the limit r s — * 0, the screened response function reduces 
to (minus) the density of states vq = mhkp / ir 2 at the Fermi 



surface. Simultaneously, c 



Up / 3, with vf = hh-p jm 



the Fermi velocity, and the Thomas-Fermi approximation for 
the dielectric function is recovered. 

Owing to its short-range interactions, the fictitious system 
can be described by Fermi liquid theoryQ. The elementary 
excitations are fermionic quasiparticles of mass m*, with an 
interaction characterized by spin symmetric (s) and spin an- 
tisymmetric (a) Landau parameters F|' a , where £ denotes the 
angular momentum channel. The Landau parameters depend 
on r s and vanish in the limit r s — ► 0, where also m* — > m. 
The speed of sound in the fictitious system can be expressed 
in these parameters as follows 



(l + F^)(m/m*)c 2 



(7) 



Due to the Pauli exclusion principle, even with an attractive 
interaction, a Fermi liquid can still support a sound mode, pro- 
vided that Fq > -1. 

Approximate calculations^^ indicate that the effective 
mass is comparable to the free electron mass, while the Lan- 
dau parameter Fq is negative. The latter implies an attractive 
quasiparticle interaction in the spin-symmetric, £ = chan- 
nel, although it derives from the Coulomb interaction, which 
repulses like charges. This unexpected finding has the same 
origin as the overscreening of a test charge. 

At the level of the first-order correction to the ground-state 
energy of a free Fermi gas, which itself is a 1-loop result, the 
origin can be understood as followsEl. Two Feynman diagrams 
contribute to the 2-loop, or Hartree-Fock correction: the di- 
rect, or Hartree term containing two fermion loops, and the 
exchange term containing only one fermion loop. Owing to 
overall charge neutrality, the direct term does not contribute 
to the ground-state energy, so that only the exchange term 
remains. Containing an odd number of fermion loops, this 
term comes with a minus sign, thus reversing the sign of the 
Coulomb interaction. Consequently, the ground-state energy 
per electron and also the pressure decrease with increasing 
coupling constant for r s < 1, eventually becoming negative. 

The inverse compressibility of the 3DEG and therefore the 
speed of sound c in the fictitious system also decrease with 
increasing r s . At a certain value r s = f s , with f s « 5.25 
according to estimatesQ, the inverse compressibility becomes 
negative and the speed of sound drops to zero, implying that 
the factor (Ar-4- Fq)/to* should vanish. The approximate 
calculations!}!! 7 ] indicate that the quasiparticle mass increases 
only slightly with increasing r$, while the values of Fq for 
r s = 2, 3, 4 show a tendency towards —1 around r s = f s . We 
conjecture that precisely at this point, Fq = — 1. In Fermi 
liquid theoryQ, this value of the Landau parameter Fq, where 
Xsc(0, 0) diverges, signifies the onset of instability, with the 
homogeneous ground state becoming unstable towards den- 
sity fluctuations. 



IV. EXPONENTIALLY DAMPED CDW 

The vanishing of the sound mode of the fictitious system af- 
fects the plasmon spectrum. Thepgeneral condition for plasma 
oscillations at a frequency lj isB e(q, uS) = 0. At zero fre- 
quency, or energy, the condition reduces to 



g 2 e(g,0) = 0, 



(8) 



where an additional factor q 2 is included for convenience. A 
physical solution q(r s ) of this condition with a positive real 
part denotes a time-independent, i.e., frozen-in modulation of 
the charge density. 

Consider condition (||) first at weak coupling, where the di- 
electric function reduces to the form (|j) in the limit of long 
wavelengths. This gives q 2 (r s ) = —ui^/c 2 , leading to a 
purely imaginary wavevector as a solution, and a screening 
length A = c/uj v \. In the limit r s — > 0, this formula reduces 
to the Thomas-Fermi result 



\/a={ira/Ar s ) 1 ' 2 , 



(4/9tt) 



1/3 



(9) 



Being proportional to the inverse square root of the coupling 
constant, the screening length (||) measured in units of the in- 
terparticle distance a becomes infinite as r s approaches zero. 

At r s = f s , where the dielectric function becomes nega- 
tive for small wavevectors, resulting in the overscreening of a 
test charge, the solution of the condition ([$]) changes qualita- 
tively. In the region r s > f s , the plasmon spectrum initially 
decreases with increasing wavevector,_uJitil reaching a (posi- 
tive) minimum after which it increases!!! 

The unique character of the point r s = f s can also be noted 
when considering the spatial average of the electrostatic po- 
tential generated by the static test charge at the origin 



d 3 Xip{x) = (f(q = 0) 



(10) 



where tp(q) = 4irZ/e(q, 0)q 2 - When the inverse compress- 
ibility becomes negative, the overall potential changes sign as 
well. 

To investigate the strong-coupling regime, we use a 
parametrized expression for the local-field correction G(q) as 
it appears in the generalized random phase approximation of 
the screened response function: 



Xsc(<7, 0) 



Xo(<Z,0) 



1 



(11) 



v(q)G(q) X o(q,0) 

proposed by Ichimaru and UtsumiEI. Here, xo(q, 0) denotes 
the response function of a free Fermi gas at zero frequency. 
The random phase approximation is recovered by setting G(q) 
to unity. Since only the zero-frequency response function is 
required to study the condition (g|), we can sidestep the dif- 
ficulties which arise when the frequency dependence is in- 
cluded in G(q) to arrive at the dynamic correction. The 
parametrized expression, which applies to the range < 
r s < 15j, incorporates Monte Carlo data on the ground-state 
energyEl as well as the ladder diagram calculation of the pair 
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distribution function at zero separation0. Importantly, the re- 
sulting dielectric function satisfies a number of exact bound- 
ary conditions and sum rules, including-.the compressibility 
sum rule (Q). As noted in the overviewEl these features and 
its simplicity makes the parametrized expression proposed in 
Ref.U convenient for applications. 

The compressibility sum rule determines the behavior of 
the local-field correction at long wavelengths. With the defi- 
nition 



limG(g) = 7o(r s )q 2 



(12) 



where q = q/kp, it follows from Eq. and a similar expres- 
sion for the noninteracting system with compressibility kq, 
that the coefficient 70 (r s ) is related to the compressibility via 



k 4a 

— = 1 lo{r s )r s 

K TT 



(13) 



The compressibility of a 3DEG can be-,extracted from the 
Monte Carlo data of Ceperly and AldeiEl thus fixing 70 (r s ). 
In particular, when the inverse compressibility changes sign it 
becomes 



r-\ 71-1 

70(js = -7- — - 

4a r= 



(14) 



Because the compressibility sum rule is satisfied, the 
parametrized expression for the dielectric function, 



e{q,0) = l-%)x S cfeO), 



(15) 



with Xsc given by Eq. (Jl 1|), becomes negative for small 
wavevectors also at r s = f a , as it should. 

With the expression (|1J) substituted and the left hand ex- 
panded in a Taylor series to order q 4 , Eq. (||) leads to the con- 
dition: 



a-o + a 2 q 2 + (HQ 4 = 0, 



(16) 



valid at long wavelengths. The quartic term is included be- 
cause the first term in the expansion, 



a-o 



(4a/n)r s 



1 - (4a/7r)7o (r s )r s ' 



(17) 



changes sign at r s = f s . The coefficients a 2 and 04 have 
no simple analytic representation, depending on the specific 
parametrization of the local-field correction G(q). They are 
best represented simply by their numerical values for each r s . 
The coefficients diverge at r s = f s . As a result of which, a 
small region just below f s is numerically inaccessible. 

The physical solution q(r s ) of the condition (|T6|), which is 
a quadratic equation in q 2 , remains purely imaginary in the 
entire regime where 1/k > 0, as it is at weak coupling (see 
Fig. [[]). The corresponding screening length A = i/q(r s ), 
being infinite at r s = 0, decreases with increasing coupling 
constant until it reaches a minimum A/a « 0.30 at r s « 3.2. 
Further increasing r s surprisingly increases A/a again until it 
becomes infinite once more at r s = f s w 5.25, where the in- 
verse compressibility changes sign and the sound mode of the 
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FIG. 1: Real (Re) and imaginary (Im) part of the wavevector q(r s ) 
solving the condition (|l6|). The grey region just below r s = f s ~ 
5.25 could not be accessed because of numerical instabilities in the 
expansion coefficients. 



fictitious system vanishes. The wavevector solving the condi- 
tion at long wavelengths is zero and the Coulomb interaction 
unscreened, as it was at r s = 0, according to this definition of 
the screening length. 

Remarkably, plotting the inverse dielectric function at this 
value of r s , indicates, to within numerical accuracy, no disper- 
sion for wavevectors smaller than the Fermi wavevector and 
l/e(q, 0) = in this range < q < hp. 

For r s > f s , the physical solution q(r s ) becomes genuinely 
complex, having an imaginary and a positive real part as well. 
The latter implies the onset of instability, where the homoge- 
neous ground state becomes unstable towards a periodic spa- 
tial modulation of the charge density. Owing to the imaginary 
part, fluctuations are, however, still exponentially screened. 
The ground state at negative compressibility is referred to in 
this paper as an exponentially damped CDW since it combines 
a homogeneous ground state with exponential screening and 
a CDW. The real part of the wavelength of the exponentially 
damped CDW is infinite at r s = f s and decreases with in- 
creasing r s in the remaining range where the parametrized ex- 
pression for e(q, 0) on which our analysis is based applies, 
i.e., r s < 15. The screening length A also starts at infinity 
and initially decreases with increasing r s . But, after reaching 
a minimum A/aw 0.54 at r s r* 8.37, it increases again. 

Apart from being exponentially damped, a CDW in a 
3DEG differs from a CDW appearing in solids also in that 
its wavevector is not necessarily given by q = 2fcp, but 
varies with r s (see Fig. [TJ). In a solid, the CDW arises due 
to electron-phonon interactions, and the charge modulation is 
accompanied.by a periodic lattice distortion both of wavevec- 
tor q = 2/spEj. 

The screening length A determined by the imaginary part 
of the wavevector solving condition ( |lq ) does not monoton- 
ically decrease with increasing coupling constant. Since the 
screening mechanism is expected to become more effective 
at stronger coupling to minimize the effect of the increas- 
ing Coulomb interaction, A fails to provide a proper mea- 
sure for this mechanism. A more relevant length scale is the 
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short-range screening length A s , in which the deviation of the 
screened from the unscreened Coulomb potential is measured 
at short distances from a test chargeci It is defined by writing 
the electrostatic potential generated by a static test charge at 
the origin [see below Eq. (flO[>] as 



< \ Z 



1 + 2 - 
IT 



dq 



1 



1 



sin(ga;) 



io yefeO) y q 
and expanding the right side in a Taylor series for small x, 

Z 



(18) 



if{x) 



with 



a 



x 



1/3 
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X 

A~ 



dq 



1 - 



e(q,0) 



(19) 



(20) 



Since l/e(q, 0) < 1, A s is always positive. The short-range 
screening length coincides with the Thomas-Fermi screening 
length in the limit r s — ► 0. As expected for a system cop- 
ing with an interaction that becomes increasingly stronger, A s 

— 1/2 

monotonically decreases (roughly as r s ) with increasing 
coupling constant in the entire range < r s < 15 where the 
parametrized expression for e(q, 0) applies. The ratio X s /a 
is unity around r s = 0.78. Whereas at r s = f s , its value is 
reduced to A s /a w 0.39. 

The exponentially damped CDW in a 3DEG arises in the 
nonperturbative, strong-coupling regime. It is worth consid- 
ering an example where a similar ground state arises at weak 
coupling, to access it in perturbation theory. Such an example 
is provided by a charged Bose gas. 



V. CHARGED BOSE GAS 

A free Bose and Fermi gas differ in that, owing to the Pauli 
exclusion principle, the latter can support a sound mode at 
zero temperature, whereas the former cannot. The single- 
particle excitation with the spectrum to = hq 2 /2m is therefore 
the only gapless mode available. By repeating the argument 
leading to the plasmon spectrum of a 3DEG at weak coupling, 
we obtain the spectrum uj 2 = uj^ + h 2 q A /Am 2 of a charged 
Bose gas in the limit r s — > 0, with the same expressions for 
the plasma frequency, uj^ = Aime 2 /m, and r s as for a 3DEG. 
The plasma spectrum of the charged Bose gas at weak cou- 
pling corresponds to the dielectric function 



lim e(q, 0) 

q->0 



-pi 



H 2 q 4 /4r 



(21) 



These formulas agree with the perturbative results first ob- 
tained in Ref£3 and Ref.Ea, respectively, using Bogoliubov's 
method. 

A solution of the condition (||) with the dielectric function 
(Elh is given by 



q(r s ) = (l + i)(mw p i/ft) 1/2 . 



(22) 



The resulting ground state is an exponentially damped CDW 
with both the wavelength and screening length expressed as 



A/a=l/(3r s ) 1 / 4 , 



(23) 



cf. the analogous expression (ffl for a 3DEG. 

As for a 3DEG, the exponentially damped CDW of a 
charged Bose gas in 3D has a negative compressibility as well. 
Specifically, the energy per particle given byc3 



E „ e 2 
= -A- 



N 2a r& 3 / 4 

to the lowest order in the loop expansion leads to 

1 5A e 2 n 



16 2a ra 3 /4 ' 



(24) 



(25) 



with A m 0.8031. A negative compressibility is possibly a 
generic characteristic of an exponentially damped CDW. 



VI. 2DEG 

Next, 2DEGs with a 1/x Coulomb potential and average 
interparticle distance a = (ir/n) 1 ' 2 are considered. The shift 
to 2D is facilitated by replacing the 3D Fourier transform 
47re 2 /<7 2 of the 1/x potential with 2ne 2 /q in the above equa- 
tions. Doing so leads to a plasma frequency 



w p i = 2irne 2 q/m, 



(26) 



which depends on q and tends to zero for vanishing wavevec- 
tors. That is, although harder than the gapless modes of the 
corresponding fictitious fermionic and bosonic systems, the 
resulting plasma modes are, unlike their 3D counterparts, gap- 
less. Most of the striking ^features of the 3D systems, such 
as negative compressibility^, the Landau parameter Fq taking 
the value — 1-at the point where thejitiverse compressibility 
changes signO, and overscreeningliio are nevertheless also 
found in 2D. 

One noteworthy difference is that the condition (|8|) in 2D 
with the weak-coupling expression for the dielectric function 
^) and the plasma frequency I 



q 2 + 2ime 2 q/mc 2 = 0, 



(27) 



has no physical solution. In accordance with the gapless 
plasmon spectrum, this implies the absence of Thomas-Fermi 
screening. 

The 2DEG's dielectric function e(q, 0) has not been de- 
termined to the extent the 3DEG's functio»-|has. However, 
the local-field correction proposed in Ref.o, incorporating 
the variational Monte Carlo data on the ground»state energy 
and the compressibility by Tanatar and Ceperleyu, serves our 
purposes as the corresponding dielectric function satisfies the 
compressibility sum rule. Specifically, with the definition 



lim G(q) = 7o(r s )g 



(28) 
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FIG. 2: Real (Re) and imaginary (Im) part of the wavevector q(r a ) 
solving the condition (f30|). 



appropriate for 2D, it follows from the compressibility sum 
rule (Q), with Xsc(q, 0) given by Eq. (|TT]), that the coefficient 
70 (r s ) is fixed by the compressibility via 



hi 



V2 



lo(r s )r s 



(29) 



The local-field correction was determined in Ref.ES only at 
discrete values of r s . We use a simple interpolating procedure 
to obtain G(q) for arbitrary values of r s in the entire interval 

< r s < 40. 

At the value r s = f s where the inverse compressibil- 
ity changes sign, with r s m 2.03 according to the Monte 
Carlo dataQ, the dielectric function becomes negative for small 
wavevectors. As for a 3DEG, the long-wavelength solution 
of the condition (0), with a factor q included instead of q 2 , 



changes here qualitatively. The resulting equation is quadratic 
in q rather than q 2 , 



a + aiq + a 2 q 



0, 



with 



a 



{V2/ir)r s 



1 - (^/7r) 7o (r s )r s 



(30) 



(31) 



while the two remaining coefficients are again best repre- 
sented by their numerical values. 

For r s < f s , no physical solution is found, implying that in 
this entire regime screening is absent and the plasmon mode 
gapless, as in the weak-coupling limit [see below Eq. (^7j)]. 
For r s > f s , a complex solution with positive real and imagi- 
nary parts emerges, signalling an exponentially damped CDW 
in a 2DEG (see Fig. g). Contrary to the 3D case, the imagi- 
nary part of the solution is larger than the real part in the entire 
regime where the parametrization applies. 

In conclusion, 2DEGs and 3DEGs at negative compressibil- 
ity, where test charges are overscreened, were argued to have 
an exponentially damped CDW as ground state. The wavevec- 
tor characterizing this state is complex and varies with the 
electron number density. The real part vanishes above the crit- 
ical density, where the inverse compressibility changes sign 
and the system becomes homogeneous. 
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